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Resummation of the chiral expansion is necessary to make accurate contact with cur- 
rent lattice simulation results of full QCD. Resummation techniques including relativistic 
CNj . formulations of chiral effective field theory and finite-range regularization (FRR) tech- 
niques are reviewed, with an emphasis on using lattice simulation results to constrain the 
parameters of the chiral expansion. We illustrate how the chiral extrapolation problem 
has been solved and use FRR techniques to identify the power-counting regime (PCR) of 
chiral perturbation theory. To fourth-order in the expansion at the 1% tolerance level, 
we find < m, < 0.18 GeV for the PCR, extending only a small distance beyond the 
^ I physical pion mass. 

<D : 

! 1. INTRODUCTION 

_> ■ 

Dynamical chiral symmetry breaking in QCD gives rise to an octet of light mesons 
recognized as the (pseudo) Goldstone bosons of the symmetry. These light mesons cou- 
ple strongly and give rise to a quark-mass dependence of hadron observables which is 
nonanalytic in the quark mass. The Adelaide Group has played a leading role in em- 
phasizing the role of this physics in the chiral extrapolation of lattice simulation results [ 
d 12 13 IH El El • Most work in the baryon sector has focused on the nucleon and Delta 
masses and the magnetic moments of octet baryons. 

The established, model-independent approach to chiral effective field theory is that 
of power counting, the foundation of chiral perturbation theory (xPT). However, this 
requires one to work in a regime of pion mass where the next term in the truncated series 
expansion makes a contribution that is negligible. Given the emphasis on determining 
the nucleon mass to 1%, such neglected contributions must be constrained to the fraction 
of a percent level. As there is no attempt to model the higher-order terms of the chiral 
expansion, one simply obtains the wrong answer if one works outside this region. 

There is now growing recognition that some form of resummation of the chiral expansion 
is necessary in order to make contact with lattice simulation results of full QCD, where the 
effects of dynamical fermions are incorporated. As we will demonstrate in the following, 
the quark masses accessible with today's algorithms and supercomputers lie well outside 
the regime of %PT in its standard form. This situation is unlikely to change significantly 
until it becomes possible to directly simulate QCD on the lattice within twice the squared 
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physical mass of the pion and with suitably large lattice volumes. Within this range, 
knowledge of the first few terms of the chiral expansion would be sufficient as higher- 
order terms of the expansion are small simply because m n is small. Straightforward 
application of the truncated expansion of xPT would then be possible. However, one 
might wonder if the lattice techniques that would allow simulations at light masses within 
2m^, might also allow a calculation directly at the physical pion mass, obliterating the 
chiral extrapolation problem altogether. Unfortunately such simulations are not possible 
in the foreseeable future, and one must resort to more advanced techniques. 

The resummation of the chiral expansion induced through the introduction of a finite- 
range cutoff in the momentum-integrals of meson-loop diagrams is perhaps the best known 
resummation method [ EJ IH EJ El III El El UD] " ~ for alternative proposals, see also Refs. [ 
ITT) IT2*] . However, there are many ways of regularizing the loop integrals which lead to a 
resummation of the expansion. For example, relativistic formulations lead to expressions 
for the nucleon self-energy which include higher-order nonanalytic terms beyond the or- 
der one is calculating [EH]- In some cases, the expressions have the desirable property of 
becoming smooth and slowly varying as the pion mass becomes moderately large [ I14j. 
Indeed all observables calculated on the lattice display a smooth slowly-varying depen- 
dence on the quark mass, presenting a clear signal that higher order terms of the chiral 
expansion sum to zero as the pion mass becomes moderately large. 

The ability to draw a curve through the lightest nucleon mass points has led some to 
argue that it is essential to adopt a relativistic approach in calculating the chiral expansion 
[ El IT5] . However, one might have some concern that a relativistic expression is so 
important in what is supposed to be a low- energy effective field theory. Indeed, a Taylor 
expansion of the relativistic expression produces exactly the same leading-nonanalytic 
(LNA) term as obtained in the heavy-baryon formulation of chiral perturbation theory 
(HB^PT). Upon including the Delta in HB^PT the next-to-leading nonanalytic (NLNA) 
structure, m^logm,,- is also obtained. 

In the relativistic formulation, an m\ correction is also generated at one-loop order. 
However, the rapid convergence of the Taylor series generated by the full relativistic loop 
calculation [ 13.J emphasizes that the relativistic corrections are small. One should not as- 
sume that the higher-order terms of the heavy-baryon expansion are included completely. 
Rather a particular resummation of the chiral expansion has been obtained which allows 
the self-energy of the nucleon to evolve slowly as the pion mass grows, thus mimicking 
the lattice simulation results. 

Perhaps the most astonishing discovery in FRR chiral effective field theory, is that the 
term-by-term details of the higher-order chiral expansion are largely irrelevant in describ- 
ing the chiral extrapolation of simulation results. Figure ^ displays the extrapolation [E] 
of lattice simulation results of full QCD from the CP-PACS collaboration [QSJ. A variety 
of finite-range regularizations are illustrated, including dipole, monopole, Gaussian and 
theta-f unction regulators, as described in greater detail below. FRR chiral effective field 
theory is mathematically equivalent to %PT to any finite order and in this case we work 
to order m\ logm^. 

The curves are indistinguishable and produce physical nucleon masses which differ by 
less than 0.1%. This is despite the fact that the coefficients of the higher-order terms 
[m\ and beyond) appearing in the FRR expressions differ significantly. For example, the 
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Figure 1. Extrapolation of CP-PACS collaboration simulation results [El to the chiral 
limit using finite-range regularization [ |3J. Differences between the illustrated dipole, 
monopole, Gaussian and theta-function regulators cannot be resolved on this scale. The 
one-standard deviation error bound for the dipole extrapolation is also illustrated. 

coefficient of the m% term vanishes for the theta-function regulator, whereas the monopole, 
dipole and Gaussian regulators provide —44, —34 and —30 GeV -4 respectively. What 
these expansions do have in common is that the loop-integrals vanish as the quark masses 
grow large. The contribution of any individual higher-order term is largely irrelevant. 
The only thing that really counts is that there are other terms that enter to ensure the 
sum of all terms of the loop integral approaches zero, in accord with what is observed 
in lattice QCD calculations. Of course, this beautiful feature of FRR expansions would 
be lost if one were to truncate the expansion at any finite order. Resummation of chiral 
effective field theory is essential to solving the chiral extrapolation problem. 

In each of the aforementioned regulators, there is a single parameter governing the range 
of the loop-momentum cutoff. It is a parameter which does not appear in standard chiral 
perturbation theory and as outlined in detail below it plays no role in the terms of the 
FRR expansion to the order one is working. However, it does play a role in determining 
the coefficients of the higher-order terms, beyond that order. As discussed in detail in 
Ref . [ E] , it governs the convergence properties of the chiral expansion. Given the level 
of agreement between the curves associated with different regulators, it is sufficient to 
parameterize the remainder of the chiral expansion in terms of this single parameter, A, 
and still maintain fraction of a percent accuracy. 

It is essential that A remain finite. In the limit A — > oo one recovers the standard 
chiral expansion of %PT truncated at the next nonanalytic term beyond the order one is 
working, spoiling the resummation of the expansion. A must be constrained to values that 



ensure that the loop-momenta contributing to loop integrals are relevant to the low-energy 
effective field theory. Indeed it is the erroneous contributions from large loop-momenta 
beyond the realm of the low-energy effective field theory that spoils the convergence 
properties of traditional [ 8J. The difficulty with this approach is that the exact 

value of A is not specified. Because one must work to finite order in the expansion, there 
is a residual dependence on A. The common approach is to vary A over reasonable values 
and include the observed variation in the final uncertainty estimate. 

2. FORMAL BACKGROUND 

In the usual formulation of effective field theory, the nucleon mass as a function of the 
pion mass (given that oc m q ) has the formal expansion: 

M N = a + a 2 m\ + a±m\ + a 6 m£ + xJ-n + XnA^A + XtadAad • (1) 

In practice, for current applications to lattice QCD, the parameters, a«, should be deter- 
mined by fitting to the lattice results themselves. The additional terms, Xnln, X^aI^a 
and XtadAad, are loop corrections involving the (Goldstone) pion. As these terms involve 
the coupling constants in the chiral limit, which are essentially model-independent [I18j. 
the only additional complication they add is that the ultra-violet behavior of the loop 
integrals must be regulated in some way. 

Traditionally, one uses dimensional regularization which (after infinite renormalization 
of cto and a 2 ) leaves only the non-analytic terms — clna"^ and cnlna""4 ^( m -n 1 1 1 ) ■> re ~ 
spectively. Within dimensional regularization one then arrives at a truncated power series 
for the chiral expansion. To fourth order 

M N = c + c 2 ml + c LNA ml + c±m\ + c N lna«4 l°g — + c 6 m l > ( 2 ) 

where the bare parameters, a*, have been replaced by the finite, renormalized coefficients, 
q. Through the chiral logarithm one has an additional mass scale, //, but the dependence 
on this is eliminated by matching C4 to "data" (in this case lattice QCD). The prime on c' 6 
denotes that this term is beyond the order to which we are working, and therefore should 
not be expected to be independent of regularization. 

In line with the implicit /^-dependence of the coefficients in the familiar dimensionally 
regulated x?T, the systematic FRR expansion of the nucleon mass is: 

M N = clq + a^ml + a% m 4 n + a$m% + xJA m ^^ A ) + XttA^a^, A) + XtadAad^, A) (3) 

where the dependence on the shape of the regulator is implicit. The dependence on the 
value of A and the choice of regulator is eliminated, to the order of the series expansion, 
by fitting the coefficients, a„, to lattice QCD data. 

To illustrate how the A-dependence of the chiral expansion is removed to the order one 
is working, we review the renormalization procedure. To leading one-loop order 

M N = a + a 2 ml + Xnln , (4) 



where x-k is the LNA coefficient of the nucleon mass expansion, 
3 
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and I w denotes the relevant loop integral. In the heavy baryon limit, this integral over 
pion momentum is given by 



/ dk yT—-i- ( 6 ) 
7r J k l + 

This integral suffers from a cubic divergence for large momentum. The infrared behavior of 
this integral gives the leading nonanalytic correction to the nucleon mass. This arises from 
the pole in the pion propagator at complex momentum k = im^ and will be determined 
independent of how the ultraviolet behavior of the integral is treated. Rearranging Eq. (jHJ) 
we see that the pole contribution can be isolated from the divergent part 

^^(f-^+ljT*^. (7) 
The final term converges and is given simply by 

where we now recognize the choice of normalization of the loop integral, defined such that 
the coefficient of the LNA term is set to unity. This choice is purely convention and allows 
for a much more transparent presentation of the differences in the chiral expansion with 
various regularization schemes. 

In the most basic form of renormalization we could simply imagine absorbing the infinite 
contributions arising from the first term in Eq. (J2J) into a redefinition of the coefficients 
ao and a 2 in Eq. (J2J). This solution is simply a minimal subtraction scheme and the 
renormalized expansion can be given without making reference to an explicit scale, 

M N = c + c 2 m 2 n + x^ml , (9) 

with the renormalized coefficients defined by 

co = CL + Xn- / dkk 2 , c 2 = a 2 - Xn- / dk . (10) 
TT Jo n Jo 

Equation therefore encodes the complete quark mass expansion of the nucleon mass 
to 0(m^). This result will be precisely equivalent to any form of minimal subtraction 
scheme, where all the ultraviolet behavior is absorbed into the two leading coefficients 
of the expansion. Such a minimal subtraction scheme is characteristic of the commonly 
implemented dimensional regularization. 

This formulation has implied that the description of the pion-cloud effects are equally 
well described at all momentum scales. As we know that chiral effective field theory is 
only valid in the low momentum regime, a minimal subtraction scheme will describe high 
momentum pion modes which have no connection with physical reality. This incorrect 
short-distance physics must be removed through the analytic terms of the expansion, 
spoiling the convergence of the expansion [ 8J. 

The pion-nucleon system in the real world is characterized by a momentum scale associ- 
ated with the size of the pion-cloud source. The axial form factor of the nucleon provides 
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a radius of {V 2 )a{L — 0.66 fm [EH], suggesting that the characteristic scale is of order 
i/( r2 )axfai ~ 0.3-0.4 GeV. Hence one should expect an improvement in the convergence 
properties of the chiral expansion upon implementing a regularization which suppresses 
the ultraviolet behavior of the loop integral. We refer to any such scheme as finite-range 
regularization (FRR). 

It should be noted that a minimal subtraction scheme does offer the advantage that 
all finite energy scales beyond the pseudo-Goldstone boson mass are integrated out. The 
finite size of the nucleon will therefore emerge as the nucleon is dressed by higher order 
pion cloud dressing. The cost of removing the explicit dependence on this physical energy 
scale is that the expansion will therefore only be reliable for pion masses below this 
characteristic scale m n < 0.3-0.4 GeV. 

We now describe the chiral expansion within finite-range regularization, where the cut- 
off scale remains explicit. In particular, we highlight the mathematical equivalence of 
FRR and dimensional regularization in the low energy regime. We introduce a functional 
cutoff, u(k), defined such that the loop integral is ultraviolet finite, 
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, k fc4 ^) , (11) 

Jo k 2 + ml 

To preserve the infrared behavior of the loop integral, the regulator is defined to be unity 
as k — > 0. For demonstrative purposes, we choose a dipole regulator u(k) = (l + k 2 /A 2 )~ 2 , 
giving 

DIP A*K + 4m7rA + A2) _ 5A 3 _ ^ m 4 m) 

16(m w + A) 4 16 i 6 m - + m - 16A 7r + '"' 1 ] 

The first few terms of the Taylor series expansion, as shown, provide the relevant renor- 
malisation of the low-energy terms. The renormalized expansion in FRR is therefore 
precisely equivalent to Eq. Q up to 0(m^) where the leading renormalized coefficients 
are given by 

A 3 5A 
c = a + Xttt^ > c 2 = a 2 - Xttt^ • (13) 
16 16 

As cio and a 2 are fit parameters, the value A takes is irrelevant and plays no role in the 
expansion to the order one is working; in this case m\. Retaining the full form of Eq. ()12|) 
will therefore build in a resummation of higher-order terms in the chiral series. 

It is straight forward to extend this procedure to next-to-leading nonanalytic order, 
explicitly including all terms up to m 2 ~ m\. Most importantly, there are nonanalytic 
contributions of order m^logm,,- arising from the A-baryon and tadpole loop contribu- 
tions. The latter arises from the expansion of the 0(m q ) chiral Lagrangian, with a coupling 
proportional to the renormalized m q expansion coefficient, c 2 . In the heavy baryon limit 

3 lQg\ f°° J; k 4 u 2 {k) 
167r 2 f 2 9 J uj(k)(A + uj(k)) 

3 2 f f 00 „ / 2k 2 u 2 (k) , , 

Xtad/tad = -T^if C2 ^U V V^+^ (15) 
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where uj(k) = ^Jk 2 + m\ and A = 292 MeV is the physical A-iV mass splitting. The 
finite-range regulator u(k) is taken to be either a sharp theta-function cut-off, a dipole, 
a monopole or finally a Gaussian. These regulators have very different shapes, with the 
only common feature being that they suppress the integrand for momenta greater than 
A. The coefficient of Eq. (JUjl reproduces the empirical width of the A resonance. 

In Eq. ()15|) . to is defined such that the term in braces vanishes at m n = 0. to is a local 
counter term introduced in FRR to ensure a linear relation for the renormalization of c^. 
This approach provides a relatively simple way to ensure that the coefficient of the tadpole 
diagram contribution is proportional to the renormalized coefficient, c%, as opposed to a 
combination of unrenormalized coefficients as done in Ref. [ I10j. Nonanalytic terms of 
the chiral expansion must have renormalized coefficients in order to recover the standard 
expansion of xPT upon taking the finite-range regulator parameter to infinity. 

The key feature of finite-range regularization is the presence of an additional adjustable 
regulator parameter which provides an opportunity to suppress short distance physics from 
the loop integrals of effective field theory. As emphasized in Eq. © by the superscripts A, 
the unrenormalized coefficients of the analytic terms of the FRR expansion are regulator- 
parameter dependent. However, the large m n behavior of the loop integrals and the 
residual expansion (the sum of the af terms) are remarkably different. Whereas the 
residual expansion will encounter a power divergence, the FRR loop integrals will tend to 
zero as a power of A/m n , as m n becomes large. Thus, A provides an opportunity to govern 
the convergence properties of the residual expansion and thus the FRR chiral expansion. 

Since hadron masses are observed to be smooth, almost linear functions of for quark 
masses near and beyond the strange quark mass, it should be possible to find values for 
the regulator-range parameter, A, such that the coefficients and higher are truly small. 
In this case the convergence properties of the residual expansion, and the loop expansion 
are excellent and their truncation benign. Note that A is not selected to approximate the 
higher order terms of the chiral expansion. These terms simply sum to zero in the region 
of large quark mass and the details of exactly how each of the terms enter the sum are 
largely irrelevant. 

The FRR expansion is summarized in Eqs. Q, ©, (fTTj). (fPIj) and (JTSJ) . Constraining 
the parameters a^ 2 4 6 to lattice simulation results provides the extrapolations of Fig. [TJ 
Table ^ summarizes the associated expansion coefficients obtained from various regulator 
fits to lattice data. Whereas uncertainties in Ref. [E] were reported at the 95% confidence 
level, we report standard errors herein. The convergence of the residual expansion as rep- 
resented by the coefficients a^ 2 4 is excellent for the FRR expansions and contrasts the 
coefficients of minimal subtraction as represented by dimensional regularization. More- 
over, the level of agreement among the renormalized quantities M^, cq and c<i for the 
FRR results is remarkable, rendering systematic errors associated with the shape of the 
regulator at the fraction of a percent level. 

In this light, it is interesting to consider the importance of two- loop contributions to the 
chiral expansion [El- The leading nonanalytic behavior is and Tr^logm^. However, 
we have already emphasized that the inclusion of contributions leads only to small 
differences in the extrapolation function, as illustrated in Fig. [T] and supported by Mjv, 
Co and C2 of Tabled Whereas the coefficient of the m% term, X5, vanishes for the theta- 
function regulator, the monopole, dipole and Gaussian regulators provide xa — —44, —34 
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Table 1 

Expansion coefficients and the corresponding extrapolated nucleon mass obtained from 
various regulator fits to lattice data. All quantities are in units of appropriate powers of 
GeV. Errors are statistical in origin arising from lattice data. Deviations in the central 
values indicate systematic errors associated with the chiral extrapolation. 





Bare Coefficients 




Renormalized Coefficients 




Regulator 


n A 
a 


4 




A 


Co 


c 2 


c 4 


M N 


Monopole 


1.74 


0.30 


-0.49 


0.5 


0.923(65) 


2.45(34) 


20.5(15) 


0.960(58) 


Dipole 


1.30 


0.37 


-0.49 


0.8 


0.922(65) 


2.49(33) 


18.9(15) 


0.959(58) 


Gaussian 


1.17 


1.22 


-0.50 


0.6 


0.923(65) 


2.48(34) 


18.3(15) 


0.960(58) 


Sharp cutoff 


1.06 


0.56 


-0.55 


0.4 


0.923(65) 


2.61(33) 


11.8(13) 


0.961(58) 


Dim. Reg. 


0.79 


4.15 


+8.92 




0.875(56) 


3.14(25) 


7.2(8) 


0.923(51) 



and —30 GeV -4 respectively. The point is that by the time the pion mass is sufficiently 
large to allow this term to contribute significantly, other complementary terms conspire 
to make the sum of terms of the FRR loop integral small. Indeed some interplay between 
C4 and X5 is observed in Table ^ where C4 displays a correlation with the magnitude of 
X5- Terms at fifth order will reveal themselves in precision effective field theory studies 
well within the power-counting regime of x?T, where m® contributions are small relative 
to the and m^logm^ terms of interest. Again, at larger pion masses, all higher-order 
terms become large and given the lattice QCD results, these higher-order terms must sum 
to zero. 

The uncertainty associated with the coefficient of has also been highlighted by 
Beane in Ref. [Ej. The traditional chiral expansion was demonstrated to show a strong 
sensitivity to the value of this coefficient. By accepting a 20% error associated with the 
truncation of the series, he suggested that the power counting regime might be restricted 
to within < 300 MeV. This is supported by Table where one can see that the 
renormalized low-energy coefficients of the expansion, Cj, have been sacrificed in applying 
the truncated expansion of to the range < m\ < 1.0 GeV 2 . 

Figure 121 further emphasizes the need to know the power counting regime of x?T. 
As one moves the regulator parameter away from the regime of 1 GeV, the otherwise 
good convergence of the residual expansion is lost, and the truncation of the expansion 
is no longer benign. The naive application of the minimal subtraction scheme outside 
the power-counting regime leads to a steeper slope in the nucleon mass relevant to the 
pi-nucleon sigma term, and produces an artificially large value. 

3. POWER-COUNTING REGIME 

Because the chiral expansion of x?T is truncated with no attempt to estimate the con- 
tribution of higher-order terms, knowing the power-counting regime of a given truncation 
is absolutely essential to extracting correct physics from %PT. One simply obtains the 
wrong answer if one works outside the power-counting regime. 

Now that the renormalized low-energy coefficients have been determined to the order we 
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Figure 2. Dipole regulator fits to lattice data for various dipole regularization scales, 
A = 0.5, 0.8, 1.1, 2.0, and oo GeV, interpolating between FRR and minimal subtraction. 



are working, we can use this information to estimate the power-counting regime of x?T. 
In the following, we quantify the power-counting regime by identifying the range of pion 
masses whereby the expansion is truly independent of regularization scheme and hence 
independent of contributions from higher-order terms in the expansion. As discussed in 
the introduction and in detail in Sec. 121 surrounding Eq. (J13)) . the FRR chiral expansion 
is mathematically equivalent to that of to the finite order one is working. In other 
words, these terms of the FRR expansion are independent of the regulator parameter, 
and any value of A is allowed. 

Fig. |21 illustrates the fourth-order chiral expansion for various dipole regulator parame- 
ters A, ranging between 0.5 GeV and oo. This therefore describes a continuous transition 
between finite-range regularization and a minimal subtraction scheme. The renormalized 
coefficients and nonanalytic terms to fourth order are automatically independent of A. 
The observed changes in the curves are simply a reflection of the changes in terms beyond 
the order one is working. Of course, only the curves with A ~ 1 GeV have a residual ex- 
pansion with good convergence properties allowing truncation of the expansion, without 
a loss of accuracy. 

The power counting regime can be identified as the regime where the curves are inde- 
pendent of A at some level of precision. As discussed earlier, this must be at the fraction 
of a percent level in order to predict the nucleon mass to an accuracy of one percent. The 
regime within which the curves agree within one percent is < m n < 0.18 GeV. This 
is displayed in Fig. 01 where the relative error between the two extremal regularization 
scales shows the sensitivity to higher orders in the chiral expansion. 
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Figure 3. With the low-energy parameters c , c 2 and c 4 fixed to those obtained by the 
fit to lattice data, the chiral expansion is shown for various values of the dipole regulator 
scale, A = 0.5, 0.8, 1.1, 2.0 and oo GeV. The left panel corresponds to the parameters q 
determined by the A = 0.8 GeV dipole regulator fit and the right panel to those determined 
in the A — > oo minimal subtraction limit. 



This sensitivity to higher-order terms also indicates that tolerance levels of 5% and 10% 
can be achieved up to pion masses of 0.25 and 0.29 GeV, respectively. These estimates of 
the chiral regime are similar to those reported by Beane [I21j. as discussed above. 

The results of Fig.0]are robust to the choice of A used to obtain the leading renormalized 
low-energy coefficients c , c 2 and c 4 of the chiral expansion. Figure El also displays results 
in the extreme (and undesirable) case where the coefficients are fixed in the A — > oo 
(minimal subtraction) limit. The associated power-counting regime is displayed in Fig. 0] 
by the dashed curve, revealing a negligible difference. 

4. CONCLUSIONS 

As we have strictly shown, FRR is mathematically equivalent to minimal subtrac- 
tion schemes such as dimensional regularization to any finite order one wishes to work. 
When working within the power-counting regime as defined above, the FRR and minimal- 
subtraction schemes are model independent. 

However, when used beyond the power-counting regime, the truncated expansion of 
%PT fails. The signature of this failure is that the low-energy coefficients of the expansion 
are sacrificed to fit the lattice results, leading to values that differ from the correct values 
of the low-energy theory. Knowledge of the extent of the power-counting regime is as 
important as knowledge of the finite-order chiral expansion itself. 

Methods to determine the power counting regime have been described herein. If one 
limits the systematic errors at the fraction of a percent level to allow a prediction of 
the nucleon mass at the 1% level, the power-counting regime is < m 2 < 0.03 GeV 2 or 
< m n < 0.18 GeV. Unfortunately, today's lattice QCD results do not come close to 
this power-counting regime of fourth-order ;^PT. With the regime extending marginally 
beyond the physical pion mass, one might wonder what role, if any, traditional will 
play in the extrapolation of lattice QCD results. The problem is exacerbated by the need 
to have sufficient data inside the power-counting regime to determine the fit parameters. 
In time, algorithms which enable dynamical- fermion simulations below m n = 0.18 GeV, 
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Figure 4. For fixed low-energy coefficients Co, C2 and C4, the relative difference in the 
nucleon mass expansion for two extremal regularization scales, where = Mm {A = 

00) — Mn (A = 0.5). The solid and dashed curve corresponds to the values of q determined 
from the fit to lattice results with dipole of scale A = 0.8 and 00 GeV, respectively. 

are likely to permit a direct simulation at m n = 0.140 GeV. 

In contrast, FRR schemes take on a model dependence when used outside the power- 
counting regime. The model-dependence associated with the shape of the regulator is at 
the fraction of a percent level, even over the range < m 2 n < 1.0 GeV 2 . As the expansion 
is truncated there is a residual dependence on the regulator parameter, A, which can 
be quantified by varying A over a range where convergence of the residual expansion 
remains acceptable. The success of FRR is realized in that this residual degree of model- 
dependence in demonstrably small, such that the chiral extrapolation of modern lattice 
results is under control. This remarkable result is extremely important in that it is likely 
to provide the only reliable method to extract physical hadron properties from lattice 
simulation results for many years. 
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